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Summary. — We present a likelihood approach to modeling multi-dimensional 
GRB _Epfc-fluence-redshift data that naturally incorporates instrument detection 
thresholds. The treatment of instrument thresholds is essential for analyzing ev- 
idence for GRB evolution. The method described here compares the data to a 
uniform jet model, in which the jet parameters are allowed to vary with redshift. 
Data from different experiments may be modeled jointly. In addition, BATSE data 
(for which no redshift information is available) may be incorporated by ascribing to 
each event a likelihood derived from the full model by integrating the probability 
density over the unknown redshift. The loss of redshift information is mitigated by 
the large number of available bursts. 

We discuss the implementation of the method, and validation of it using simulated 
data. 

PACS 98.70.Rz - gamma-ray sources; gamma-ray bursts. 



1. Introduction 

Analyses of HETE-2- and BeppoSAX-localized Gamma-ray bursts with known red- 
shifts provide evidence that the isotropic-equivalent GRB energy Z?i SO evolves with red- 
shift 1 . The assessment of the strength of this evidence is linked with the issue of 
detection thresholds, which can create apparent energy-redshift correlations by trunca- 
tion. The modeling of such thresholds is a complex and delicate matter 

We present here a method, based on the likelihood function, for analyzing GRB E p k~ 
fluence-redshift data in the presence of instrument thresholds. We have developed the 
method with a view to addressing the question of whether GRB properties do, in fact, 
evolve with redshift. 

2. The Method 

We compare observations of GRB data with a model of GRB source properties, using 
the likelihood function. 

The data consists of GRB fluences and E p k values, with errors, and GRB afterglow 
redshifts, where available. We use the results of fits of the Band function to integrated 
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GRB spectra. The fit parameters that we use as data are y p k = \og w (E p k) and ])a = 
log 10 (^4), where A is the amplitude — the scale factor that multiplies the Band function. 
These parameters are supplemented by their errors, o^and a a, respectively. 

The model is a generalization of the variable-angle uniform jet model of Lamb, Don- 
aghy and Graziani 2 . The generalization is that the parameters of the distribution of 
jet opening angles are allowed to evolve with redshift. 

The observables of the model are: 

• 2, the redshift. 

• x iso = log 10 (£; iso /10 52 erg). 

• x pk = ^°&io(^ P k A keV), where E^^is the E p k parameter in the source frame. 

• Xj et = log 10 (r2 :(et /47r) = log 10 (l — cos9j et ), where 8j e t is the opening half-angle of 
the uniform jet. 

• x-y = log 10 (i? 7 /10 52 erg), where E 1 is the prompt gamma-ray energy actually 
emitted: x 7 = Xi SO + Xj et . 

The model is introduced through the differential event rate: 

pi^Z , X^ , X j e t 5 ^dz dx^y dxj e t dXpfo — 

(1) (j,(z)dz x Lo(x 7 )dx 7 x K(x je t\z)dx je t x C( x p k \x iso )dx^ 

where 

• is proportional to the Rowan-Robinson 0] SFR function, multiplied by (1 + 
z) dV(z)/dz, where V(z) is the proper volume of the sphere bounded by redshift 
z. 

• lu(x 7 ) = exp (^—^^-j^—^j , with T = 0.33 2 and x 1 a fit parameter. 

• n{x jet \z) = Cljet~ 5 x fijet < %et < ^fjf ^ and ^fJt ma y be functions 
of redshift. 

• C( x pk C) \ x iBo) = exp f-^-^— ^!f^!\ w ith b = 1.95 and A = 0.13, is the 
broadened Amati relation |2] ■ 

The likelihood function is given by 

^instruments Ni ( ff.ihi) (M) r (M)\l 

( 2 ) L = n n{ — r T 

;=i i=i I 1 J 

where Ni is the expected number of events observed by the Z-th instrument, given its 
specific threshold, and f(y^' % \ Dp£\ z^ 1 '^) is the individual likelihood of the i-th event 
observed by the l-th instrument. 
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In Eq. |(!2J), the expected number of events N is given by 
(3) N t = J dzdx jet dx 1 dx^ k rc) p(z, x 7 , x jet , x^^) x (% et /47r) x r}i(z, x iso , x p s k rc) ). 
Note that: 

• rji(z, Xi SO , z^pk^) * s ^ ne detection probability by the l-th instrument of a GRB whose 

vital statistics are z, Xi SOl x pk . It is here that the instrumental threshold is in- 
corporated. r\ is given by an analytic function that depends on the instrument 
thresholds. We adopt the threshold values of Band 0. 

• The factor f^ et /47r is the probability that our line-of-sight lies inside the jet. 

• Ni is reducible to a three-dimensional integral over z, Xi SO , and x^ c \ This integral 
must be performed numerically. 

In Eq. (J2J), the event likelihood / is expressed in terms of variables xa = log 10 (A) , 

(srt 
pk 



x pk = \og w (E p k), which are related to x iso , x^ c \y 



x pk C) = x p k+log w (l + z) 

fd L (z) 2 \ 

x iso = x A + (a + 2)x p k + log 10 I ^ - - J - A, 

where a is the Band-function low-energy power-law index (assumed equal to -1 through- 
out), A is a unit-conversion constant, and d^{z) is the luminosity distance. We then 
have 

f(jjA: Vpki z) — J dXj e t dXpk dx A p[z, X^ — Xiso "Ejeti %jeti %pk ) 

1 (y A - x A ) 2 1 (y pk - x pk ) 2 



x exp 



This integral may be performed analytically by standard techniques. 

In the case of events for which no redshift is known, the event likelihood should be re- 
placed in the full likelihood by its marginalized version, i.e. f(jjA, Vpk, z) — > / dz f(yA,Vpk, z). 

3. Implementation and Validation 

GRB evolution is currently implemented by allowing the three parameters of K,(xj e t\z) 
(power-law index, lower- and upper-bounds of power-law domain) to be themselves 
power-laws in 1+z. Evolution may be turned off by setting the indices of these power-laws 
to zero. 

The resulting six parameters are supplemented by a seventh — cc 7 . 

The likelihood function is interpreted according to the principles of Baycsian inference. 
It is converted to a posterior density by multiplication by an assumed prior density on 
the model parameters. The shape of the resulting probability distribution in parameter 
space may be used to estimate (or set constraints upon) the parameter values. The peak 
of the distribution constitutes the best-fit parameter estimate, while confidence regions 
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Fig. 1. - Estimating the power-law index of the Qj e t distribution. 



may be obtained by considering level curves of the posterior density containing prescribed 
amounts of probability (90%, say). 

We are currently validating the code by running it on samples of simulated (non- 
evolving) GRBs produced by the simulator described in 0. Figure ^shows the result 
of estimating the power-law index of the f2 jet distribution using 25 events drawn from a 
simulation that assumes an index 5 — 2. The true value of the index lies within the 90% 
confidence interval. 

Work on estimating x-y and the other power-law parameters (which are more broadly 
distributed and more highly correlated) is in progress. We will be applying the experience 
obtained from analyzing simulations to real data from HETE, BeppoSAX and BATSE 
in the near future. 
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